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Abstract 

The nonlinear integral equations describing the spectra of the left and right (con¬ 
tinuous) quantum KdV equations on the cylinder are derived from integrable 
lattice field theories, which turn out to allow the Bethe Ansatz equations of a 
twisted “spin —1/2” chain. A very useful mapping to the more common non¬ 
linear integral equation of the twisted continuous spin +1/2 chain is found. The 
diagonalization of the transfer matrix is performed. The vacua sector is analysed 
in detail detecting the primary states of the minimal conformal models and giv¬ 
ing integral expressions for the eigenvalues of the transfer matrix. Contact with 
the seminal papers ram by Bazhanov, Lukyanov and Zamolodchikov is realised. 
General expressions for the eigenvalues of the infinite-dimensional abelian algebra 
of local integrals of motion are given and explicitly calculated at the free fermion 
point. 
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1 Introduction 


An integrable system in 1 + 1 dimensions is characterized by the existence of a number 
of integrals of motion in involution equal to the number of degrees of freedom. Con¬ 
sequently, the exact computation (of the eigenvalues) of these commuting observables 
is a very important problem. In this respect, an important and widely used method 
goes through (an integrable discretisation of the model and) the diagonalization of the 
discrete integrals of motion by a Bethe Ansatz technique. At the end one is left with the 
problem of finding the continuous limit of the eigenvalues of the (discretised) theory. In 
order to achieve this goal, a powerful analytic method was introduced by Destri and de 
Vega for the periodic sine-Gordon theory on a cylinder J5[. Although we will follow the 
way depicted in [3], we want to mention that similar results were first obtained by pfj 
in the case of the critical behaviour of the 6- and 19-vertex models by using a different 
technique. We do prefer to use the procedure of Destri and de Vega since it has been 
successfully applied to the calculation of ground state energy pj and first excited states 
energies jB] of sine-Gordon model, being our aim in the second paper the treatment 
of massive perturbation of Conformal Field Theories (cfr. j7] as review paper on the 
subject). 

In the paper |B] we have formulated a generalisation of the Algebraic Bethe Ansatz 
to diagonalise the transfer matrices of two lattice models, obtained after discretisation 
(on a space lattice) of two - left and right - quantum (m)KdV equations. The lattice 
formulation of one lattice (m)KdV was already performed by Kundu P(, stressing the 
necessity of modifying the usual Yang-Baxter relation to take into account the non¬ 
ultralocality of the theory. The integrability structure was found in a braided version 
of the Yang-Baxter relation, which has been also attracting mathematical interest ra 
in itself. Afterwards, we have found the Bethe equations, which may be thought of as 
those of a “spin —1/2” (“wrong” sign) XXZ chain with a twist. The origin of the 
twist in the framework of |5j is in some sense dynamical , i.e. automatically generated 
by the theory and depending on the number of Bethe roots, almost like in m and 
unlike in [|B]. It is also worth mentioning that a space and time discretisation of the 
KdV was proposed in [El and analysed by using the Baxter TQ —system approach Ex 
these discrete-discrete integrable systems are indeed important also for their continuous 
parents [ Til . 

In this article we address the problem of writing nonlinear integral equations gov¬ 
erning the spectrum of the left and right (continuous) quantum (m)KdV theories. In 
other words, the eigenvalues of the transfer matrix - which encodes all the local and 
nonlocal conserved charges - are expressed as nonlinear functionals of the solutions of 
the aforementioned nonlinear integral equations. Besides the importance in Statistical 
Field Theory, they can be also interpreted as explicit expressions for spectral determi¬ 
nants of VT - symmetric spectral problems [T31ITB] . This procedure will eventually give 
the eigenvalues of the local integrals of motion of the quantum KdV theory. After El, 
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it became well known that these local integrals of motion represent one series of local 
commuting charges of conformal held theories without extended symmetries m ms 
(cfr. [EHHHI] for details concerning the matrix Lax formulation). In addition, this se¬ 
ries is exactly that which is still conserved, after suitable modification, if the conformal 
held theory is perturbed by its $( 1 , 3 ) operator. Actually, applying the investigations of 
the present paper and of [ 8 ] to the setup elaborated in [[TQj about the KdV theory, it 
is only a matter of calculations to extend the subsequent results to the only other series 
ESI, characterising the only other integrable perturbations. The very intriguing feature 
of the conformal formulation in terms of integrable hierarchies is that the conformal 
monodromy matrix already contains the perturbing held and the screening operator 
nsi- A similar discretisation should be also allowed in conformal held theories with 
extended symmetry algebra (e.g. for W algebras it is of basic importance the setting 
of [2T]; cfr. {22 for development in a peculiar case). 

We remark that similar results have been obtained in pioneering works mia by 
Bazhanov, Lukyanov and Zamolodchikov from a different starting point and via a dif¬ 
ferent approach. They define directly in the continuous held theory the transfer matrix 
as an operator series acting on Virasoro modules (for 0 < /3 2 < 1/2). On the other 
hand, we intend to create a bridge between conformal held theories and the powerful 
algebraic Bethe Ansatz formulation of lattice KdV - based on this interesting 

generalisation of the Yang-Baxter algebra —, showing that the continuum held theory 
limit can be recovered in the braided case as well. As particular consequence, in the 
second paper we will be able to formulate the $( 1 , 3 ) perturbation of conformal held 
theory within the same framework of [3J by merging two KdV theories. In any case, 
the monodromy matrix of dei is different from the naive continuous limit jSJ of our 
monodromy matrices, which, unlike in DEI. are not solutions of the usual Yang-Baxter 
equation. 

In section 2 we summarise the main results on (m)KdV theories obtained in JH]. 
In section 3, starting from the Bethe equations of i.8], we write the nonlinear integral 
equation for each vacuum of the left and right quantum KdV equations. In section 4 we 
use this equation to calculate up to quadratures the energy of the vacua. In section 5 we 
apply the nonlinear integral equation to find exact expressions for the continuous limit 
of the vacuum eigenvalues of the transfer matrix. From the asymptotic expansion of 
these eigenvalues we obtain the local abelian charges, mentioned before as characterising 
the < h(i, 3 ) perturbation. In addition, we have solved exactly the theory in two cases: at 
the free fermion point /3 2 = | and in the limit of infinite twist (a) —> + 00 ). Some results 
are also compared with results and conjectures in DEt In Section 6 we summarise our 
work describing also possible further applications. 
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2 Quantum lattice (m)KdV theory 


In a previous paper [Hi we considered the left and right mKdV equation, respectively 

,, 3 o . 1 /// „ 3 n , 1 ttt . « 

d T v = -vv + -v , d f v = -vv + -v , ( 2 . 1 ) 

for the fields v = —p' and v = —p', dehned as spatial derivatives of quasi-periodic 
Darboux helds p , (p in the interval y, y G [0, R]. The quantizations of the Darboux 
fields are the Feigin-Fuks bosons 0, 0 (23j . which satisfy the commutation relations 

= ( 2 . 2 ) 

where 0 < 0 2 < 1 and s(z) is the quasi-periodic extension of the sign function: 

s(z) = 2n + 1 , n < z < n + 1 ; s(n) = 2 n , n G Z . 

In terms of a discretisation of the Feigin-Fuks bosons, 0 m = 0 , 0 m = 0 ( m 2 ^), 

we dehned the ^-periodic operators V± living on a IV-site lattice with spacing A and 
length R = NA: 


V m = 2 (02m-l — 02777+1) + (02)71-2 


i) — (02m—1 — 02)77+l) + (02)71-2 — 02m) (2.3) 


V + = - 

m 2 |_ 


( 02 m —1 — 02 m + l ) + ( 02 m -2 — 02 m ) — ( 02 m -1 — 02 m + l ) + ( 02 m .-2 — 02 m ) -(2-4) 


These operators are the quantum counterparts of the discretisation of the mKdV vari¬ 
ables u and h EU) and satisfy the non-ultralocal commutation relations, first introduced 
in 0 , 


[K+, K + ] = 
\ y „ 7 ki = 

WZ 7 V~] = 


iTT/3 2 

2 

m/3 2 

2~ 

m(5 2 


{& m — l,n 3 m , n — l ) ) 

(fim— 1,77 *^ 777 , 77 — 1 ) 7 

( fim — 1,77 2 S mn T $ m , n — l ) 


(2.5) 


In terms of we now dehne (adopting a different normalisation with respect to 
| 8 ]) the /e/t and right conformal monodromy matrices 


M(a) = e 2 L N (a)... Li{a) , (2.6) 

_ X7T _ _ 

M(a) = e — L N (a) ... L 1 (a ), (2.7) 

where left Lax operators L m (a ) and right Lax operators L m (a) are given by 


/ e Ae"e* y? " 

L m (a) = e 4 I ^ a e -iv+ e iv- 


,L m {ot) = e 4 




Ae “e 


Ae 


( 2 . 8 ) 
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As a consequence of non-ultralocal relations (1231) . monodromy matrices mm 
satisfy braided Yang-Baxter relations (see, for instance, relation (4.4) of |SJ), first intro¬ 
duced in [2] and ra- The commuting transfer matrices associated to them generate, 
restoring the continuum, one series of local commuting charges of, respectively, left and 
right conformal field theories. 

Despite the complication coming from the presence of a braided Yang-Baxter rela¬ 
tion, we have derived in |S] the Bethe equations for the monodromy matrices QUID- 
They read: 


sinh(o: s — a r + m(3 2 ) 
JLi sinh(ci' s - a r - i-K/3 2 ) 

r^s 


sinli (a s — — £ In 

sinli (a s + yy — e In 


N/2 


e ei7r/3 2 (-f+2i) . 


(2.9) 


The equations may be formally interpreted as being those of a twisted “spin —1/2” 
XXZ chain, for a sign is the only difference with respect to the usual twisted spin +1/2 
chain equations. Indeed, l is the number of the Bethe roots and the parameter e is 
equal to —1 for the left conformal case (monodromy matrix (12.61) f and to +1 for the 
right conformal case (monodromy matrix (12.71) 1. The connection with notations used 
in formulae (6.35, 6.47) of H is given by q = e~ m ^ 2 and A r = e OLr . 

In the same paper fHj we wrote also the eigenvalues of the left and right lattice trans¬ 
fer matrices on the Bethe states. Using the same notations as in m such eigenvalues 
read 

A at (a) = A%(a) + A^(a), ( 2 . 10 ) 


where 


and 


A^(a) = e 


— 0 =F£*7T (3 2 l 


1 sinh(a: — ay ± iirf3 2 ) ± 


n 

r —1 


sinli (a — ay) 


Pn(<*) 


p N (a) = e 


inJ3^N_ (‘IR 

— o I 11 4 


Vlv e 


„,+ 2 . [ inf] 2 R' 

1 smh I ea ± £— -In — 


N/2 


( 2 , 11 ) 


( 2 . 12 ) 


In this paper we will perform the continuous limit of these quantum lattice (m)KdV 
theories, with the aim of studying quantum (m)KdV theories. 


3 The nonlinear integral equation for the vacua sec¬ 
tor 


We want to write a nonlinear integral equation equivalent to Bethe equations m f ° r 
the particular solution minimising the energy (the vacuum). We introduce the analytic 
function in the strip |Ima;| < min {</, 7r — C} ? 0<C <7r (l 11 is taken in the fundamental 
branch) 




sinh(i£ + x) 
sinli — x ) ’ 


(3.1) 
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to define, when N G 4N, the counting function, 

Z N (x ) = y 0 (x-e In y, <f>( x - a r , vr/3 2 ) + 27 t/? 2 £ ^ + y ^ , (3.2) 

in terms of which the Bethe equations m have the logarithmic form 

Z N (a s ) = tt(2 I a -l- 1), I S EZ. (3.3) 


First of all, we are interested in the vacuum solutions of the Bethe equations: all the a s 
are real and Z^{a s ) are equal to all the numbers of the form m between Z^{— oo) 
and Zn{+ oo), i.e. there are no holes (cfr. sub-section 5.2 for further details). It follows 
from the definition (EH that these two extremes are given by 

Zjv(±Oo) = ±y (it - 7 T/3 2 ) ± l(7 T - 27 T/3 2 ) + 2n/3 2 £ (l + y^ . (3.4) 

This solution minimises, for the range of values of (5 2 

| < /J 2 < |, (3.5) 

the left (e = —1) and right (e = 1) hamiltonian eigenvalues, defined as logarithmic 
derivatives of the transfer matrix eigenvalues, respectively 


H e ({a r },R,N) = 


T 


Ae 


2 sin 7 t(3 2 


—2 ea. 


9 r-( 

95 Ga,( “ 


i R ~ f~ 

a —£ In N 2 


—2eot 


d_ 

da 


G+(a) 


a= £ ln 


(3.6) 


where 

G%(a) = ln[p^(a) _1 Aj\r(a)]. (3.7) 

On the other hand, the vacuum solution is known to minimise in the region (13.51) 
the lattice sine-Gordon Hamiltonian eigenvalues H lsg({ov} , R, N) [21]; left and right 
hamiltonians eigenvalues (ED can be regarded as the conformal limit of H: 


HJ{a r },R,N) = lim H LSG ({a r - £0 - £ In R/N}, R, N). (3.8) 

©^+oo 

Supposing /3 2 in the region ED, we now want to write a nonlinear integral equation 
for the vacuum counting function in the continuous limit: 


N —> Too , R fixed . 


(3.9) 


Besides, this limit means, because of ED, that also l 
put 


fd 2 = 


P_ 

pi 


Too and therefore, even if we 


(3.10) 
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with p < p' coprimes, it turns out to be ill-defined on the exponential factor of Q. 
Therefore, we slightly modify the method in m and we need to consider different ways 
for l to go to oo according to the sub-sequences 

l = np + k , 0 < k < pJ — 1, (3-11) 


where n —> +oo and n is kept fixed. The utility of this limit procedure is that we can 
redefine the counting function m up to suitable multiples of 2tt to obtain 

Zn(x) = (x - eln-^, ( x ~ “n tt/3 2 ) +2vt^£ (j5 2 - l) +2tt£lo , (3.12) 

where u is finite in the continuous limit (the double braces denote the fractional part): 


u = 



(3.13) 


The presence of the constant divergent term in definition (ETT2D will be made clear in 
the following. 

Now, we can use standard techniques (see 0, 0, [H|) to rewrite a sum over real 
Bethe roots of a function with no poles on the real axis as follows 1 : 


/jv(ov) 

r =1 



/ “boo rj nr r • 

—/^(x)Imln 1 + e iZ *r( x + i0 ) 

-oo 27T L 


(3.14) 


The use of (HUD in mu gives 


(x) = y^{ x ~ £ ln +27 F t £ (^ 2 ~~ X ) + 27T£uj + 


/ 2tt^( X ~ Vl 7T ^ Zn ^ ~ 2 / _ V' -/? 2 ) Im ln 


2tt 


(3.15) 

]_ _j_ e iZ N (y+iQ) 


After defining: 
1 


K(x) = —</>' (x, vr^ 2 ) , $(x) = </> ^x, ^/? 2 ^) , L n (x) = Im ln 


1 + e 


iZjv(ai-t-iO) 


(3.16) 


the Fourier transform of relation mrsi) is 

A N $(k)e~ i£kln % n K(k) f ... ^ 9 

Zjv(A:) = —---2-—Tjv(A:) + 47t 2 £ 

2 1 -K(k) 1 -K(k) 

where we have used the notation: 


u; + 


N 


(y- 1 ) 


m 


r+oo 


f(k) = / dxe~‘ kx f(x) 


J 1 - A'(fc) ’ 

(3.17) 


(3.18) 


1 We are fixing l to be even in the following (i.e. e lZ N( a s ) _ _iy as the case l odd can be treated 
along the same lines and in the end does not add new information. 
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Explicitly, the Fourier transforms of the first two functions in (18.1611 


are: 


$(fc) = - 


2 vr 


sinli k | 


ik 

With the definitions 
7 = tt(1-/3 2 ), G(k) = 


sinli jk 


K(k ) 


, £(*) = 


sinli ^ (f — 7r/? 2 ) 


sinli 


(3.19) 


sinli k — 7 ) 


1 — /\’(A') 2 cosh /c| sinli k f — 7 ) 

equation (HUD can be written in the following way: 


, 07=^07, (3.20) 


N &(k)e~ i£kln % 

Z N (k) = — - y h 2G(k)L N (k) + 2 tt£ 

^ 1 — A (/c) 


~ N L 1 
a, + TM 1 -y 


S(k). (3.21) 


We now want to Fourier antitransform and to perform the continuous limit. There¬ 
fore, it is necessary to work out in such a limit the behaviour of the TV-dependent 
term 


F n (x) = j 


_ f + °° dk , kx / N$(k)e 


—iek In 


2tt 


ttN r +0 ° dk 


2i J -00 2nk 


2 1 - K(k) 

>»«). sink* (f - f/J 2 ) 


sinli k^/3 2 cosh k — ^fl 2 ^ 


(3.22) 


This expression obviously diverges in the continuous limit. To define a finite count¬ 
ing function in such a limit we perform a TV-dependent shift Aat in x and define the 
conformal counting function Z(x) as 


Z(x) = lim Z N (x + A N ). 

N —►+00 
NA=R 


(3.23) 


We have to choose Aat in such a way that 


F N (x + A n ) 


ttTV r+°° dk cik ( x _ eln R +An ) sinli k (f - f/? 2 ) 

2i 7-oo 2itk sinh/c|/ 3 2 coshA; — |/ 3 2 ) 


(3.24) 


is finite in the continuous limit. A convenient way to evaluate (18.2411 is the residues 
method. The poles of the integrand lie at k = 0 and on the imaginary axis. Since we 
are in the region (13.51) . the poles on the imaginary axis with smaller modulus are at 
k = ± 771 / 7 . We choose A N as follows 


R 7 

- e In — + A n = x- In 
N 7r 


A . 7 xf3 2 

N 8111 2(1 — /3 2 ) 


(3.25) 


where x can be equal to ±1 and A is an arbitrary positive constant which reflects the 
translation invariance in x of conformal field theories. 






















If x — — 1, then in the continuous limit — £ In + An —> +oo: hence, to evaluate the 
integral for fixed x we can choose a contour closing in the upper ^-complex half-plane 
and containing a small semicircle surrounding the pole at k — 0. Since the contribution 
from the imaginary poles are damped by the exponential factor containing An, the 
leading term is given by the residue at the pole with smallest modulus, i.e. k = in/^, 
minus the contribution of the semicircle around k = 0: 

N —*• +oo , F n (x + A n ) — — TVe^V 11 ^ + — lj + o(N°) , (3.26) 

where o(z ) means “order less than z”. 

If x — I; then in the continuous limit — elnjl + An —> — oo and we consider a 
contour closing in the lower /e-complex half-plane and containing a small semicircle 
surrounding the pole at k — 0. As before, we take into account only the contribution 
from the residue at the pole with smallest modulus, k = —itt/x, and we subtract the 
contribution from the semicircle around k = 0: 

N— >+oo, F N {x + A N ) = Ne~ x e hl £ - - l'j + o(N°). (3.27) 

Collecting these formulas we get 

N —> +oo , F n (x + A n ) = A X e^ x - X -j7r - i) + o(N °). (3.28) 

We now make the choice 

X = -e, (3-29) 

and insert (13.2811 in the continuous limit of the Fourier antitransform of (13.211) shifted by 
An- The constant divergent term disappears and we obtain for the conformal counting 
function (13.231) the following equation: 


Z(x ) = — Aee 


+ 2 


/•-boo 


dy G(x — 2 /)Im In 1 + e 


o iZ(y+i0) 


+ £UJ , 


(3.30) 


where G(x) is the Fourier antitrasform of G{k ) (13.201) 


G{x) 


r+ °°dk Jkx sinli fc (f — 7 ) 


2vr 


2 cosh k | sinh k 


(3.31) 


Remark 1 In writing (13.171) we have paid the necessary attention to the zero modes. 
Besides, we could more correctly derive equation (THUD working out our calculations 
on the derivative of Z(x) |25j . but we think that this derivation is simpler, once all the 
subtleties are under control. 

Remark 2 Although equation (13.301) has been found for 1/3 < ( 3 2 < 2/3, yet it can 
be considered without problems in the whole region 0 < (3 2 < 1, defining in such way, 
by analytic continuation, a state in the continuous limit which minimises the energy. 
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Remark 3 In fact, equation (13.301) coincides with the nonlinear integral equation 
for the spin +1/2 XXZ chain which shows Bethe equations (12.91) with the l.h.s. inversed 
and (3 2 = Hence, the mapping (3 2 —> 1 — /3 2 relates the twisted spin +1/2 and spin 
— 1/2 XXZ chains. This connection between spin +1/2 and spin —1/2 XXZ chains was 
already conjectured in m, but only in the case of infinite volume. 

Remark 4 When e = —1, we can make contact with Bazhanov, Lukyanov and 
Zamolodchikov’s paper j2p Indeed, equation (13.301) coincides with equation (3.8) of J2j, 
derived for 0 < (3 2 < 1/2, if we choose A = 2M cos > where 


M 



T(1 -+)W^ , 


(3.32) 


and if we identify p of t 2| with ujj 2. Such an A is positive, as it has to be, in the whole 
region 0 < f3 2 < 1. Therefore, in spite of the differences highlighted in the Introduction, 
(13.301) with this choice of A can be regarded as an extension of (3.8) of to the domain 
0 < + < 1 . 

Equation (ECU) is the nonlinear integral equation describing the cu-vacuum of the 
continuous theory. Its counting function solution allows us to evaluate the infinite 
continuous limit of a sum over the Bethe roots of the c+vacuum resetting formula 
(13.141) as 


i 


lim V/W(<+•) = 

N —>+oo ^ 

NA=R r=1 


f+°° dx 


r+°° d x 


' —oo 27T 

where we have defined the continuous limit 


—f'(x)Z(x) + 2 / —f(x )Imln 1 + e iZ(x+i0) 

-oo 27T L 


, (3.33) 


f(x)= lim f N (x + A n ). 

N —>- + oo 

NA=R 


(3.34) 


Relation (13.331) contains two different contributions; we will call them, for the sake 
of brevity, the bulk term and the finite size correction term. They can be identified 
inserting in the first term of the r.li.s of (13.331) the expression for Z(x) coming from 

(EHH1) : 


1 r + OO rjx 

lim J2f N (®r) = 2 —f'(x)Imhi 

K-1 J ~°° 277 


1 + e 


iZ(x+iO) 


r+°° dx 

-oo 27t“ 


r+oo 


f{x)\Aee £ ~< x — 2 / dyG(x — y )Imln 


+ 

X _|_ e iZ (y+ i0 ) 


(3.35) 


— ELU 


We can arrange this expression 


^ r+oo Ax 

iim X)/jv(ar) = fb + 2 / ^/+) 


r+oo 


N —>+oo 
NA=R r=1 


dy[S(x — y)—G(x — y )]Imln 1 + e iZ (y +l0 ') 5 

(3.36) 
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where fb indicates the bulk term 


/ 4 -OO sjnp 


Aee 


ecu 


(3.37) 


and the other term gives usually the finite size correction. It is useful to write these 
terms using the Fourier transform of /, defined in EUBf 


lim X!M a r) = fb 


N —►+00 . 

NA=R r=1 


r+o o d jQ 

/ —In 

/ —oo 7r 


l + e 


iZ(x-\-iO) 


where 


J f (x) = [ + °° —e ikx ikJ(k)f(k ), J{k) = -- 

A 7-00 2vr 1 h y ’ 2 cosh k^ sinh & (f “ 2 ) 


(3.38) 


(3.39) 


In the next two sections we will evaluate the bulk term and the finite size correction 
term of the conformal hamiltonian and of the eigenvalues of the transfer matrix, both 
in the left (chiral) and in the right (antichiral) case. 


4 Finite size correction to the energy of the vacuum 

The aim of this section is to extract, from the finite size correction to the vacuum 
energy, the conformal charge of the theory and the conformal weight, characteristic of 
the cb-vacuum. 

Let us start by writing the energy (EH) of a Bethe state: 


H £ ({a r },R,N) = J2 h e(a r ,R,N), 


r =1 


where 

h e (x,R,N) = 


sA e 


3 — inf } 2 A — e)—x 


-,— in 0 1 R — e)—x 


2A sin vr (3 2 \ [sinh (4 - e In § + sinh (x - e In f + 

+ 


+ 


4 2 --4 2 ]}. (4.1) 

In the expression - see (ETBgjl . (ECU - for the hnite size correction of the energy what 
enters is the Fourier transform of EH. Explicitly we have, using 


(4.2) 


„ g-ifce In ^ gin 7rfc ^ 2 ( -jik 

ikh £ (k , R, N) = 2ei xk— -—-sinh nkfl 2 - eind 2 

v ’ ’ ; A sin 7 x(3 2 sin A ^ a 2 M 

Therefore, the quantity Jh{x) (13.3911 in this case is given by: 

J h {x) = lim 


f +oa dk ik • 

x -( l -/3 2 ) e\n 

A ■ tt/3 2 1 

iv sm 57T^ 

1 Nenk sinh 

Tikf3 2 — — emfi 2 ^ 

/-oo 2vr e 


i?sin7r/3 2 cosh k | 
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This integral can be evaluated using the residue method. Since we are going to the 
limit N —> Too, we close the contour of integration in the upper /^-complex half-plane 
if £ = 1 and in the lower fc-complex half-plane if £ = — 1. In both cases the only 
non-vanishing contribution in the continuous limit is given by the pole with smallest 
modulus, i.e. that at k = ieir/'j, and hence turns out to be 


Jh(x) = 


2/1 sin ^§ 2 ) sin (ef ■ggr) 
R (1 — /3 2 ) 2 sin 7T/3 2 


— £ — X 

e 'y 


(4.3) 


Let us indicate by E^(R) the finite size correction to the energy of the u- vacuum, 
defined as the second term of the r.li.s. of (|3.38j) with Jh given by ED- The effective 
central charge is proportional to the A-indipendent quantity 

6 RE q (R) _ 12A sin 2 ^ sin (ef r+oc dx _ £ g 

7T 7T (1 — /3 2 ) 2 sill TT/3 2 J- oo 71 

As well explored 0I30E1, the integral in ED can be computed using a derivative 
lemma , yielding the effective central charge 


lTe 


iZ(x-\-iO ) 


(4.4) 


c ‘ff 



(4.5) 


Now we can identify the (conformal held) theory, i.e. the central charge c(/3 2 ) indepen¬ 
dent on the twist a), and the conformal weight of the Bethe state h reinterpreting ED 


where 


and 


Ceff = C - 24h , 


c =13-6[P 2 + - 


P 2 


, p 2 u> 2 c - 1 CO 2 c - 1 

h = + —— = —- + 


(4.6) 

(4.7) 


4tt 2 24 4/3 2 ' 24 

Nevertheless, if we want to implement constraint (ism oo can take only the values: 


(4.8) 


iv — — , 0 < s < p' — 1 

pi 


(4.9) 


Therefore, inserting ED in ED we get the effective central charge for fixed P 2 = p/p': 

/ 6s 2 \ 


Ce ff ~ 1 


pp' 


(4.10) 


In the particular case of minimal unitary models, p' = p + 1, the effective central charge 
reads 

/ 6s 2 ^ 

Ce// V p(p + 1 ). 


, 0 < s < p. 


(4.11) 
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These values of c e // identify all primary states of minimal unitary models which are on 
the diagonal of the Kac table. Sublimits different from BUD may yield other primary 
states and the entire table should be completed by holes and/or complex Bethe’s roots 
excitations [28) . 


5 Continuous limit of the eigenvalues of the transfer 
matrix 


It is surely interesting to use the results of section 3 in order to evaluate the continuous 
limit of the vacuum eigenvalue of the transfer matrix 42HJD, in particular as in this 
limit the integrals of motion it generates in the asymptotic expansion are connected to 
the local conserved charges of the aforementioned conformal field theory Hi- 
Referring for notations to formulae (12 .1 0112.121) . we define 

F£(a) = lnA^(cii) , F^(a) = In A^(a ), (5.1) 


and we consider first F + . From vrnurm we have 

1 sinh(a — a r + inB 2 ) . ^ 

- -r- reiTp 

smh(a — a r ) 

in/3 2 R\ ii t/3 2 


F+(«) = 


X> 


r=1 


+ 


N 

+ 2 


In sinh \ea + £ 


In — 

N, 


, 2 R 
+ In — 

2 N 


ea 


(5.2) 


We want to evaluate for the vacuum 


F + (a) = lim Ft(a + \ N ). (5.3) 

N —> + oo 

NA=R 

The last addendum in (15.21) gives a contribution whose continuous limit is zero if 0 < 
/ 3 2 < 1/2 and infinity if 1/2 < f3 2 < 1. However, one can get rid of this divergence by 
defining F + (a) for 1/2 < /3 2 < 1 as the analytic continuation of the same function for 
0 < /3 2 < 1/2. 

With this regularisation procedure in mind, we will calculate, for 0 < /3 2 < 1/2, 


i 


F + (a ) = lim 

N—> + oo . 

NA=R l ' r '= 1 


E 


In 


sinh(a! + An — a r + z7t/ 3 2 ) 
sinh (a + A^v — cx r ) 


— iSTT/3 2 


(5.4) 


and we will release the constraint on /3 2 . As already pointed out, the constant term in 
(fPl) has to be treated indipendently. Its contribution is —i£ir/3 2 l and, since l is even, 
we can restrict ourselves to the subsequences in which the integers n and k, contained 
in BUD , are both even. With this restriction, using the relations (13.11)1) . (13.111) . the 
constant term in B3> can be written as follows 

— ien—n + 2niZ. (5.5) 

p/ 
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The multiples of 2ni can be discarded since the transfer matrix depends only on 
exp[E + ]. Remembering the definition of the twist (Id. 181b contribution (15.51) becomes 
(the double square brackets denote the integer part): 


— IE7TCU — IE7T 


pk 

pt 


The integer part gives a factor to exp [ F + ] which can be equal to ±1, depending on the 
choice of k. For the sake of simplicity we will carry on the calculations supposing that 
k is chosen in such a way that this factor is 1: the vacuum eigenvalue of the transfer 
matrix for the other values of n differs form the result we will obtain in this section 
by a global factor —1 and hence its treatment is redundant. Therefore, without losing 
generality, we can conclude that the constant term in El gives a contribution to F + 
which is equal, up to multiples of 27ri, to —ienuj. 

Using this result and (I8.88|) for dealing with El, we obtain 


F + (a) = /+(«) + / - 


+°° dx r f+°° dk 


7r L J— oo 


— e ikx ikf + (k, a)J(k) 
2tc 


Im In 


1 + e 


iZ(x-\-iO) 


isncu . 


where / fe + is a bulk contribution, J{k) is given by (18.891) and 

. , sinh (a - x + mf3 2 ) 

f (x, a) = In- —7 -s-• 

smh(a — x) 

After defining 

y = a 

we will study F + (y ) = F + (a ) in the two regions: 

REGION 1 
REGION 2 


• 7 
2 ’ 


7 7 

- < Irny < - , 
2 2 


7 7 

-7r < Imy < -. 

2 2 


(5.6) 

(5.7) 

(5.8) 

(5.9) 
(5.10) 


The behaviour of F + (a) in the other regions of the complex plane can be obtained 
using the periodicity property, 


F + (a + in) = F + (a ), 


(5.11) 


which comes from expression (15.41) for F + (a). 

Let us start with region 1. We remark that, when |Iim/| < 


f + (x, a) = - i<f> (x , (5.12) 

where (j) is defined in El- From (15.61) it follows that we have to calculate the Fourier 
transform of the derivative of f + (x,a ) as function of x. From the definition of K(x) - 
first of the (ET5D - and from its Fourier transform - second of the EH - we obtain: 


ikf + (k, a) 


-2nie~ iky 


sinh k — 
sinh k% 



(5.13) 
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This result implies that, in this case, Jf+(x,a) (13.8911 bears the form (after some sim¬ 
plifications) 

J '^-^ = hSZ ikek{ ^^ < 5 - 14 > 

This integral is finite if Irny < ^ and the result is m 

i 


Jf+{x-y ) = - 


7T 


i'f cosh - {x — y) 


(5.15) 


Now, after inserting in (15.61) . we have 

, . r+°° dx 1 

F + (a) = /+(«) + / 


i'f cosh -{x — y) 


Im In 


1 + e 


iZ(x-\-iO) 


IS7TCJ . 


(5.16) 


and we need to calculate (E3ZD 

A + («) = 


r+°° dx df + (x, a) ■ 


Aee E 'y x — euj 


d 27r dx 

Using dSH we arrive, after some manipulations, to the integral 
f ’+ 00 dx sinh inf3 2 


fb ( a ) = ~ 


J-oo 27T C osh [y — x — cosh (y — x + 
which is simple to be explicitly calculated when 0 < /3 2 < 1/2: 

7T 2 


Aee 


't — EUJ 


(5.17) 


(5.18) 


/ b + (a) = is A ctg — e £ -y v + ie'nuj , 

2y 


(5.19) 


and, instead, diverges if 1/2 < f3 2 < 1. However, in the regularisation procedure we 
chose, we define the limit (15.31) for 1/2 < [3 2 < 1 as the analytic continuation of its 
value for 0 < j3 2 < 1/2. Hence, collecting formulae (15.161) and (15.191) . we get 


7 r 


F + (a ) = ieA ctg — e e i y + 

2y 


r+°° dx 


iy cosh -{x — y) 


Im In 


1 + e 


iZ(x-\-iO) 


(5.20) 


for the whole interval 0 < /3 2 < 1 and if Irny < 

Let us now focus on region 2. The evaluation of ikf + (k,a ) can be done starting 
from m and using 


ikf+(k,a ) = 


(5.21) 


It follows that Jf+(x,a ) is now given by the expression 


Jf+(x, a) = 


r+oo 


2 J-c 


dke ik ( x ~ y -^) 


sinh k7} 


coshfc^ sinh k 


(5.22) 
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which, as far as we know, cannot be expressed through simple operations on elementary 
functions. On the other hand, the evaluation of the bulk term, defined by formula Em 
can be performed explicitly. With the help of [2^ we obtain the following result, valid 
in region 2 and again for 0 < (3 2 < 1/2: 


/ b + (a) = ieA 


3 e ”( !/+, 2) + i 


ISTIU) — IEUJ , 


sin 


2 7 


(5.23) 


With the regularisation procedure we discussed before, we can extend this result to the 
whole range of values 0 < (3 2 < 1. Putting together EH , EH with (15.61) . we finally 
obtain the expression for the value of F + (a) in the entire interval 0 < /3 2 < 1 and in 
the region 2, ^ — n < Irny < — 


F + (a) = ieA ■ 


Sill tt 
27 


3 ~e*(y+i f) 


+ / + °° — [ + °° dke ik ( x - y ~^) 


ieuj + 


i sinh k ^ 


(5.24) 


7T 7—oo 


2 cosh /),'! sinh k — 


fin In 


1 + e 


iZ(x+iO) 


This result can be also obtained as analytical continuation of the function E2DD beyond 
the singularity of the integral kernel at Imy = —7/2: in the x-plane the pole crosses the 
real axis and we need to take it into account through the residue theorem and hence 
using equation Em) . 

Let us now consider F~. From (I2TTD . (IZT2D and EU> we have that 


+ 


f n(®) = S[ ln 

r= 1 

N " 


sinh(o( — a r — m (3 2 ) 


sinh(Q: — ay) 


lnsinh (ea — A-A- — In — 

V 2 N . 


+ ien/3 z 


+ 


, 2 R m(3 2 

+ ln--ea + e— 


(5.25) 


We want to evaluate (for the vacuum) 


F (a) = lim F N (a + An). (5.26) 

IV^ + oo 

NA=R 

Again, we have that the continuous limit of the last addendum in (15.251) is zero if 
0 < (3 2 < 1/2 and infinity if 1/2 < /3 2 < 1. However, according to our regularisation 
procedure, we define F~ for 1/2 < /3 2 < 1 as the analytic continuation of its value for 
0 < (3 2 < 1/2. Hence, we have 


F~(a) 


lim 

N —» + oo 

NA=R 


E 


In 


sinh (a + Atv — ay — Z7r/? 2 ) 
sinh(a + Ajv — ay) 


+ ein/3 



Comparing E2B with El, we notice that 

F~{a ) = — F + [a — infi 2 '') , 


(5.27) 


(5.28) 
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and this means that the behaviour of F (a) in all the complex plane can be obtained 
from results (15 .201) , dE2H) and (EHTD for F + (a). We stress that, because of the period¬ 
icity property «nu). relation dE2HD can be alternatively written as 

F~(a) = —F + ( a + iy) . (5.29) 

To summarise, we have shown that the vacuum eigenvalue of the transfer matrix, 

A(a) = exp[F + (a)] + exp[— F + (a + iy)], (5.30) 

can be obtained from relations (Em (ictk and from the periodicity property (I^TTD - 


5.1 Asymptotic behaviour 


We now want to extract information on the vacuum eigenvalues of the local conserved 
charges of conformal held theories. In the classical case (/ 3 2 —> 0) these appear JDj as 
coefficients of the asymptotic expansion of the logarithm of the transfer matrix A(a) 
around the complex A = e a = oo. Likewise to !Z9 we expect that the local integrals of 
motion are encoded in F + (a). But now there are two different asymptotic expansions 
of F + (a) following from (15.2011 and (15.241) respectively. More explicitly, since 


Rey 


-eoo 


cosh -(x 
7 v 


y) 


+oo 

= 2 £( 

n=0 


■iy 




where the symbol = represents asymptotic expansion, from flEZDD we have that, when 
Rey —> — eoo in the strip Imy < 

(5.31) 


7r 


F + (a ) = ieA ctg—e £ A + 
2y 


2 + °° 
*7 “ 0 


J2(-l) n e e 5 y(2n+1) 


r+oo 


dxe 


-e^x(2n+l) 


Im In 


X e iZ(x+iO) 


whereas, from (JOl, it follows that, when Rey —> —eoo in the strip ^ — 7r < Imy < — 
F + (a ) = — ieu + ieA ——^ e~ £ ^^ y+l ^ + 

- 1 )"" 


sm 75- 
27 


(5.32) 


+00 f 

^E 


7r 2 (2n+l) 


n=0 ^ 7 COS 27 

1 ') n hn ^(n+l) 

J-J tan—— _2*_/ 


e £ y(y+if )(2n+ 1 ) f + °° dxe -^n+l) lmln h + ^+*0) 
J — OO 


+ 


e ^(v+^ )("+i) f + °° dx e~ E ^ x{n+1) lm In [l + e iZ(x+l0) 

J — OO *- 


7T — 7 


In fact, likewise to [20], our final aim is to expand asymptotically the logarithm of the 
transfer matrix 


lnA(a) = In {exp[F + («)] + exp[— F + (a + iy)]j , (5.33) 
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for y G C and 



Rey —> —eoo , 

(5.34) 

with the parameter (3 2 restricted to the semiclassical interval 



0 < /3 2 < \ ■ 

(5.35) 

We will carry on the asymptotic analysis of In A (a) dividing the 

one-period strip 

V 

2 

: — 7r < Imy < 3 2 

(5.36) 

in the following domains: 



DOMAIN 1 

: ^(3 2 - 7T < Imy < | - tt , 

(5.37) 

DOMAIN 2 

7 7 

: - - 7T < Imy < -- , 

(5.38) 

DOMAIN 3 

7 T 3 

: -y < Imy < tt-- 7 , 

(5.39) 

DOMAIN 4 

: 7T — ^7 < Imy < 3 2 . 

(5.40) 


Let us start from domain 1. The expansion of F + (a) is given by (15.3111 with y —> 
y + in, thanks to the periodicity (I5TTD . and similarly the asymptotics of — F + (a + ry) 
is equal to (15.311) with a global minus sign and with y —> y + fy. This shift compensates 
the global minus sign, and therefore we have 

2 

F + {a ) = isActg—e~ £ ^ y+in) + (5.41) 

2y 


n +oo 

_ (—l) n e £ y (y+in)(2n+l) 

1 7 n = 0 


r+oo 


dx e 


-e?-x(2n+l) 


Im In 


1 + e 


iZ(x+iO) 


— F + (a + i'y) 


. 7T —f—v 

ieA ctg— e + 

2y 

o +oo 

b„= 0 


(5.42) 


+oo 


dxe“ e ^ (2n+1) Imln 


X giZR-MO) 


Which of these two expansions dominates in the asymptotics of the logarithm of the 
transfer matrix (15.3311 is given by the comparison of the asymptotic values of their 
leading terms. Hence, 

lim Re 

Key—>—eoo 

implies that (EUD dominates: 


ieA ctg 


7r 

27 


—e—y 

e ^ 


— IE A ctg 


7T 

27 


g-ey 


= +oo 


lnA(a) = — F + (a + i , y) = 

7T 2 2 + °° 

27 


IeA ctg—e £ ~< v + — y^(—l) ra e £ i y ( 2n+1 '> 


r*+oo 


dxe 


-e^x(2n+l) 


Im In 


*7 .^o 


(5.43) 

X _)_ e iZ(x+iO) 
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In the domain 2, when Key —» —eoo the asymptotic expansion of F + (a) is given 
by formula whereas the asymptotic expansion of — F + (a + iy) is still given by 

relation (15.4211 . The leading terms of (15.4211 and (15.4211 satisfy the limit 


lim Re 

Rey—>—£00 


7r 


ieA ctg — 

2y 


—e—y 

e 


ieA ■ 


Sill tR 
27 


S) 


= +c>o, 


therefore, when Key —> —eoo, In A(a) = — F + (a + iy), i.e. in domain 2 the expansion 
of the logarithm of the transfer matrix is still given by formula (15.4411 . 


We now analyse domain 4. In this domain the asymptotic expansion of F + (a) is 
given by expression dim whereas the asymptotic expansion of — F + (a + iy) is given 
by minus expression ( 1021 ) in which y is substituted by y + iy — in. Hence, up to the 
leading term 

2 

F + (a) ~ ieA ctg — e~ £ ^ y (5.44) 


and, on the other hand, 


— F + (a + iy) ~ — ieA 


* 7T^ 

sin 7T- 




2 7 


(5.45) 


Again, one can show that in the limit Key —eoo the leading term ■ ion . is dominat¬ 
ing: 


lim Re 

Rey—eoo 


• ^ 71 
ieA ctg — 
2y 


e e ~i v + ieA 




sin 


2 7 


= +oo; 


therefore, when Key —> —eoo, lnA(a) = F + (a ) (I5.41jl and, since (15.411) coincides with 
(1021) and (1031) . the asymptotic expansion of In A(a) in the domain 4 is still given by 

(1031) . 


We are now left with domain 4. Again, the asymptotics of F + (oi) is given by 
expression (1QT1) . and, as a consequence, its leading term by (1031) . On the other 
hand, — F + (a + iy) is given by minus expression (15.411) in which y is substituted by 
y + i'y — in and, consequently, its leading term 


— F + (a + iy) ~ —ieA ctg — e £ -y( y+n 


(5.46) 


is subdominant, i.e. 


lim Re 

Rey—eoo 


70 


ieA ctg — e 
2y 


+ ieA ctg 


7r 


2y 


e ~^(v+i 7-in) 


= +oo 


Therefore, when Key —> —eoo, the asymptotic expansion of lnA(o:) is still given by 

(1031) . 
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The validity of such a picture can be partially extended to the interval 

\ < 0 2 < 1 ■ (5.47) 

Indeed, following the same procedure as in the semiclassical case, one can show that in 
the strip 

3 T 7 

--7 < Imy < - 

the asymptotic expansion of In A(a) is given by (15. 4, ‘II) : however, in the strip 

7 t 3 

- - 7T < Imy < --7 , 

which completes the in periodicity, the asymptotic development of In A (a) does not 
follow in a simple way from CT E3Z). However, it is still possible to write such an 
expression, which is found to be different from 4031)- 

To summarise, we have shown that in the limit Rey —> — eoo and in the domain 
|/3 2 — 7T < Imy < |/ 3 2 the asymptotic expansion of the logarithm of the vacuum 
eigenvalue of the transfer matrix in the semiclassical domain 0 < /3 2 < 1/2 is given by 


2 + 2 ? 

In A(a) = - ism e~^ y + - £ {-l) n e £ ^ 2n+1) c n+l I 2n+l , (5.48) 

m 71=0 


where 


7T 


m = — A ctg — , 
27 


(5.49) 


^n+1 


and 


1 + £ f3 2n+2 7T2 F ( 

(n + 


(i + O) 

(n +1)! 2 r 

(i + l 

( n + l) 

'«) 


[ T(l — / g 2 )]- (2n + 1) ( 1 +0 , f = 


P 2 


1^, 1 = [ + °° dxe~ e 5 x{2n+1) Imln 


Cn+l 7 " / -°° 


1 + e 


iZ(x+iO) 


1-P 2 ’ 

(5.50) 

(5.51) 


Of course, the asymptotic behaviour of lnA(a) in all the complex plane follows from 
(15.481) and the periodicity property (15.111) . Moreover, relation (15.481) is still valid in the 
domain 1/2 < (5 2 < 1, provided we restrict it to the strip — §7 < Imy < |. Without 
this restriction, we would have that the expansion is governed also by the dual non-local 
charges, appearing in the second series of (15.321) . 


In the next subsection we are going to prove that the local integrals of motion hn+i 
in (EUBD are also normalised by the c n as usually in Conformal Field Theories. For the 
entire interval 0 < /3 2 < 1 they are given by the asymptotic expansion (15.311) of F + (a). 
The case n — 0 


+OO 


C\ / y J —00 


dxe 


'Ini In 


l + e : 


iZ(x-\-iO) 
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can be evaluated using the derivative lemma im but an explicit formula for the cases 
n > 1 is still missing. 

Moreover, we want to stress that, when £ = — 1, formulae (l5.48H5.51j) coincide with 
Conjecture (2.21) of [2J for the regime 0 < /3 2 < 1/2. In fact, the two asymptotic 
expansions for In A (a) and their domains of validity are the same if we compare the 
quantum Lax operator (EED and the classical one in |8] to obtain the connecting A 2 = 
e 2 " = e 2y+ * 7 , and if we notice that (15.4911 reproduces formula (4.15) of |2j, 

p d 1-2/3 2 \ 

m = 2v/? V A:g( r(l^/3 2 )T», (5.52) 

\ 2 — 2 / 3 2 ) 

after choosing A as in Remark 4. On the other hand, we have found two different 
asymptotic expansions for F + (a), i.e. formula (15.3111 . valid for —7/2 < Irny < 7/2, 
which coincides with (4.14) of j2j, and formula (15.5211 . valid for 7/2 — n< Irny < —7/2. 


5.2 Limit L) — > +OO 


This limit linearises equation (TOUll into an inhomogeneous integral Wiener-Hopf equa¬ 
tion of second kind (which is solvable using standard Fourier transform techniques). I11 
fact, after introducing new normalisations 


e = -x, z(o) = z(^o), G(6) = ^g(-6 

7 \7T / 7T \7T 


(5.53) 


equation (EH) reads 


Z{9) = -Aee~ £d + 2 


r+oo 


d9'G(9 - 9')lm\n 


1 + e 


iZ(Q'+iO) 


+ £UJ . 


(5.54) 


Let us confine ourselves to the case e — —1 (the case e — 1 can be treated along 
the same lines). We will show in the following that equation (15.5411 is properly valid 
in the range 0 < Co < 7t/ 2/3 2 - although the results of Section 3 hold for all Co > 0 
thanks to their analyticity. Therefore, we calculate the exact solution of ()5.54jl in the 
limit Co — > +cx), meaning by that the analytic continuation of the solution valid in 
0 < Co < it/2(3 2 . Although the periodic solution do exist, we prove that the analytic 
solution is not periodic, for we need to modify the (15.5411 itself by imposing analytic 
continuation. Indeed, its derivation is based on the assumption that the only points 
6 on the real axis for which = —1 are the Bethe roots and this is true only 

if 0 < Co < 7t/ 2/3 2 : in this case, the holes 9^ h \ such that e l Z( eW ) — _l an d do not 
satisfy the Bethe equations (EH), all lie on the axis Im 9^ h ‘ = 7r. We see this because 
these holes are also the zeroes of the transfer matrix (lOHl and we also understand that 
when Co increases from 0 to 7 t/2 /3 2 the real part of the first hole from left decreases; at 
Co = tt/2/3 2 , the real part of that hole reaches — 00: 

lim A(0 + iO) = 0 => + e~ i02 “ = 0 ; (5.55) 

Re 6^—00 
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then for u) > 7 t/2 (3 2 this hole jumps onto the real axis and migrates to the right 1 . 
We have borrowed our reasoning from a slightly different context (2S| • Nevertheless we 
can go further in our analysis. Indeed, it is known that the presence of holes on the 
real axis for d) > 7 t/2 /? 2 implies that the integral term of .54^1 - which should count 
only contributions coming from the poles of Im In 1 + e lZ (. e '+ l °) corresponding to Bethe 
roots - in fact receives contribution also from the holes on the real axis. Therefore, in 
the presence of Nh holes on the real axis at the positions 9 1 ,... 9n h , we need to subtract 
the holes contribution as follows 6 . EB] • 


+oo 


Z(0) = Ae e + 2 / dO' G{6 - 0')hn In 


l + e : 


iZ(0'-\-i0) 


N h 


— u — i In S(9 — 9h) ■ (5.56) 


h =1 


In ()5.56j) we have introduced the sine-Gordon soliton-soliton scattering amplitude S(9), 
whose explicit expression is 


S(6) = exp 


2 J- 


+°° sin k9 sinh - h r) 

00 k cos h Fl sinh ( 


kiv 2 

2 7 


)J 


(5.57) 


and whose connection with G(6) is 

G(9) = 


i d 
2 tt d6 


In 5(0). 


Likewise, in the presence of holes the quantity F + {6) (15.20115.24(1 is modified into 


F +(h \e ) = F + (6 ) +2 iJ2 arctane ( ^- 0) , 9 = 


N h 


7r 


h =1 


7 


-V' 


(5.58) 


As a result we can hnd the exact values of uj at which one additional hole migrates onto 
the real axis. Indeed, the condition 

lim A l ' h \9 + iO) = 0 , 

Re 6 —>—00 

applied to the expression of the vacuum eigenvalue of the transfer matrix in the presence 
of holes on the real axis, 

?+Q 1 ), 


A {h \9) = exp[F +W (0)] + exp [~F +{h \9 + m )\, 


(5.59) 


is still equivalent to 

+ e= 0 , 

because in the limit Re 6 * —> — 00 the added sum in (j5.58l) produces in (15.5911 an overall 
factor e mNh . Therefore, we conclude that when ^ < & < (n + |^ ^ we have 

that n + 1 holes are moving right on the real axis. 

1 We are really indebted to R. Tateo who has suggested us the possibility of the presence of holes 

on the real axis, although our treatment would not change if the holes were not present. Yet their real 
presence is conceptually important. 


22 
















We are now going to write the nonlinear integral equation (15.561) and the eigenvalues 
(ETK31) in an alternative form which is suitable for the study of the limit Co —> +oo: the 
integration path can be shifted in order to avoid contributions from the incoming holes. 
In this respect, we make the hypothesis that there is a separator s(Co) on the real axis 
such that for all the Bethe roots 9 ^ and for all the holes on the real axis 9^ 


6 {h) < s(Co) < 9 {r) . 


(5.60) 


Under this hypothesis we can move down the integration path (—oo + iO, s(to) + iO) in 


1(9) = 2 



9') Im In 


1 +e 


iZ(0'-\-i0) 


(5.61) 


of (15.541) . in absence of holes (0 < Co < 7 T/2/? 2 ), in such a way that, when the holes reach 
the real axis, they will never cross the integration path: 


rs(u>) 


1(9) = -i / d9'G(9 - O') In 1 + e iZ(M) 


r+oo 


— I 


M £>) 


d9'G(9 - 9') In 


1 + e' 


iZ(0'+i0) 


p+oo 


+ i 


d9’G(9 - 9') In 


(5.62) 

e -iZ(0’-iO) 


This modified expression contains only contributions from the Bethe roots and holds 
by construction for any Co. Clearly, when 0 < u < 7 t/2 /3 2 it coincides with (15.611) . 
Algebraic manipulations allow us to rewrite (15.621) in the compact form 


r s(u) 


r+oo 


1(9) = / d9’G(9 - 9')Z(9') + 2 / d9'G(9 - 9 ')Im In 




1 + e 


iZ^Q'+iO) 


(5.63) 


where the first integral contains Z(9') thanks to the absence of holes and roots. Eventu¬ 
ally, using relation (15.631) we are able to extend analytically the equation (|5.54l) Vw>0: 


Z(9) 


Ae e — oo + 

/•-boo 


r s(u) 


d9' G(9 — 9')Z(9') + 


’s[u) 


d9'G(9-9') Imln 


1 + e ; 


iZ(0'+i0 ) 


(5.64) 


This form is very suitable for the study of the limit Co +cxd because all the nonlinearity 
is confined in the last term. In fact, the last term can be neglected in the limit Co —> +oo, 
for it can be approximated replacing Z(9) with Ae e — Co: 


p+OO 


I nl (9)~ 2 / d9'G(9-9 ')Imln 


’s(ui) 


1 + e 


i(Ae( 6 '+ i0 '> -6j) 


After a simple change of variable we have 

G(0- In 

L V 


I nl (9) ~ 2 / + °° dt- 

4-ii+Ae s (“) 


~A~ 


t + 00 


Im In 


1 + e 


*t(l+iO) 


(5.65) 


(5.66) 
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Now, for 9 < 8' the integral EQ) furnishes us an asymptotic expansion [8 — 6' —> — oo) 
on residues 

G[8-d') = Y J Cie ai{e ~ e, \ di> 0 . 

i 

Using this expansion we easily argue that we can neglect the last term of (15.641) in the 
limit u) —> +oo provided for instance the inferior integration limit of (15.661) 

— d> + Ae —■> +oo , as u) —> +oo . (5.67) 

This is realised by the asymptotic behaviour - A can be always chosen > 1 - 


oj —> Tex) , s(u))~lnu;, 


(5.68) 


which will be proved a posteriori - with the existence proof of a separator - using the 
explicit solution. 

In conclusion, we can approximate in the limit u —> +oo equation (15.641) with a 
linear equation of Wiener-Hopf type, 

Z[8) = Ae e -u+ f H dd' G[9 - 8')Z[8'). (5.69) 

J —OO 


We can solve this equation with the usual factorisation technique [213 finding the Fourier 
transform of 


obviously defined as 


Z-(0) = 

Z-(k) 


m 

if 6 < s[lo) 

0 

otherwise 

rs(u>) 


/ 

d6e~ ike Z (9) 


(5.70) 

(5.71) 


This integral is well defined when Irn/c > 0 because, when 8 —> —oo, Z[8) goes to a 
non-zero constant. 

With the choice of A as in Remark 4 the solution of (15.691) can be expressed by 


Z-(k) 


2 iTrz 7 r p — 

~V^1 r (| - f) r (l ~ikt) 



Inserting in fTTul i the expression of M (HO^l) . we simplify the solution: 


(5.72) 


Z-(k) 


i n ik+ § 2 ik 

n 

u 

-«¥) 

i 

k f3 2lk ~ 2 Co lk ~ 1 Y | 

(i- 

2 ) 

i r i 

4- 

ikfj 


[r(i - ,3 2 )]® <1+£ ). 


(5.73) 


We remark that it is now possible to verify a posteriori the existence of a separator and 
assumption Indeed the evaluation of the Fourier antitransform of ra . 


Z.[8) 


+oo+ir rllp . 

?-ri kd Z_(k), T> 0, 

-oo+ir Z7T 


(5.74) 
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can be approximately performed, in the limit Co —> +oo, applying the residues technique 
after closing the contour of integration in the lower /c-complex half-plane. The leading 
contributions come from the poles at k = 0 and at k = —2 i/(1 + £): 


Z_(0) ~ —2[i 2 tb + 2 jt 3/2 -+!—+ (- 

r(++)A 


‘00 


2/3—1 


= 2(l-/3 2 )0 


(5.75) 


Now the condition Z_(9) = (2 n + l)7r identifies the holes, which happen to appear at 
values of 9 


0 ^ 


In 


r(| + /3 2 ) (paV- 2 ? 


2 - 2(p 2ir 3 / 2 r(l - IP) \ 2ir j 

Hence, it follows that, when Co —> +oo, all these holes go like: 

9 n ~ In u). 

Therefore property dESHl) holds for the separator s(u;). 


(2u7T + 7T + 2(3 2 oo) 


(5.76) 


We now want to consider the vacuum eigenvalues of the transfer matrix and in 
particular the conserved charges defined by (ED) , eventually finding the exact form 
(15.5011 for the normalisation coefficients c n+ 1 . Formula (15.511) reads 


(‘n-i -1 f 271-i- I 


r+oo 


d9e e(2n+l) Imln 


1 +e 




and it is valid when 0 < uo < ti/2/3 2 . With the described procedure, we can write the 
expression valid also for Co > ti/2/3 2 , i.e. when holes are present on the real axis: 


Cn++ 2 n+l 


1 /■«(&) 
2 J—oo 


/* + OO 


dde { - 2n+1)e Z(0)+ / d9 e (2n+1)6l Im In 
J s(uj) 


1 + e' 


iZ(O-hiO) 


(5.77) 


As well as for the nonlinear integral equation, in the limit oo —+ +oo the second adden¬ 
dum in the r.li.s. of (EZZP is negligible, yielding 


i.e. from dSUD 

And finally yields 

^71+1-^272+1 — 


c n+l I 2n+1 ~\ [ SH d9e^ e Z(9), 

Z J — OO 


Cn+ihn+l - -Z_(i(2n + 1)). 


(5.78) 


(5.79) 


i + e f3 4n+A u 2n+2 r( 

(n + 

0 

(i + O) 

(n + 1)! 7r 2ri+ 5 2 2n + 3 r 

(! + ' 

( n + i) 

'«) 


[T(l-+)]” (2n+1)(1+ + (5.80) 
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On the other hand, local conserved charges of conformal held theories can be constructed 
in terms of Virasoro generators |18j . Their vacuum eigenvalues can be calculated for 
the first cases and we easily understand that they behave for large u) as 


hn+\ - 



2n+2 


(5.81) 


The normalisation coefficients c n+ 1 of expansion (EHKl) are independent on uj. Therefore, 
their exact value is obtained by inserting in <EBZ9 the expression (j5.81j) for 2n+1 • 
Therefore we obtain, in conclusion, the anticipated result (|53Up . 


5.3 The case (3 2 = 1/2 

When f3 2 = 1/2 (free fermion point) the nonlinear integral equation (18.3011 becomes 
trivial, because - as follows from definition (|3 .3 111 - G(x ) =0. Its solution is 

Z(x) = -Aee~ 2£X + e2-kuj . (5.82) 

In this case it is possible to evaluate explicitly the eigenvalue of the transfer matrix. 
We fix the arbitrary constant A as in Remark 4 (but we do not choose e) and we first 
consider the case in which the spectral parameter a is in the region 1 (Q: 

0 < Ima < 7 t/2 . 

Some care is necessary, for the bulk terms of formulae ( EH ) and diverge logarith¬ 

mically when (3 2 = 1/2. The analytic regularisation is again very suitable and as shown 
for instance in ED it is equivalent to a minimal renormalisation scheme: we define 
j3 2 = 1/2 — ?7, and we will let 77 — 0, maintaining only the finite part 2 . Considering 
first F + (a ), we have from (15.2011 

F+(a) = B + (a) + C + (a), (5.83) 


where the bulk divergent part reads 


B + (a) 


—b 27 t( 1 + "0(1) — hr 7r) + 47 xea 

V 


e -2ea + } 


(5.84) 


where 0(1) = — C = —0.57721566490... (C is Euler-Mascheroni constant) is the renor¬ 
malisation constant depending on the regularisation scheme (for instance, if we had 
considered the dependence of the function (15.8411 on y , we would have obtained a dif¬ 
ferent constant). Instead, the regular finite size correction is given by 


j /*+ 00 p 2ecx. 

C + (a ) = — dk— -- 

7T Jo k 2 — e~ 4ea L 


ln(l + e 


2tt 2 ie 2lQ k—2 txiuj\ 


— ln(l + e" 


—27 v 2 ie 2l0 /c+27mj5 


(5.85) 


2 More in general, divergences of the bulk terms occur for (3 2 = 1 — l/2n, with n > 1. 
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We now go to the limit 7] —» 0 using a minimal renormalisation scheme: 


B + (a ) = [27t(1 + -0(1) — In it) + Att ea] e 2sa . 

On the contrary C + {a) (15.851) is finite and can be explicitly calculated 

C + (a) = — slnT ^ + u + £7re~ 2EOl '\ + elnT + uj — £Tre~ 2ea 


(5.86) 


einu — [27t(1 — In 7r) + Att£oi — i£7T 2 ]e 


2] _ —2ea 


(5.87) 


Therefore, gathering the previous two formulas, we obtain, in the region i (EHD, 

F + (a ) = B + (a ) + C + (a ) = — elnT + to + £7re _2£Q ^ + (5.88) 

+ £ In T ^ + uj — £7re~ 2£Q ^ + 27T0(l)e~ 2£ “ + 77r 2 £e _2£Q — £77ru;. 

For what concerns F _ (a), it comes from the general relation (15.281) . At P 2 = 1/2, when 
a varies in region 1, a — ii\P 2 varies exactly in region 2 (—7r/2 < Ima < 0). Therefore, 
the value of F“(a) is the finite part of the ij —> 0 limit of minus expression (lA2dl) in 
which a is shifted by —iirP 2 . Performing this limit as we did before for F + , we obtain 
the renormalised values 



F (a) = B (a) + C (a) + 2£ittu) , 

(5.89) 

where 

B~ (a) = B + (a) - 2TT 2 i£e~ 2sa , 

(5.90) 

and 

C~{a) = C + (a). 

(5.91) 

As a consequence, 



f-(«) = 

—£ In r Q + a; + £tt e~ 2ec ^j + £ In T Q + a; - £7re _2£ “) + 


+ 

27T0(l)e _2£ " — i'K 2 £e~ 2ea + £77ru;. 

(5.92) 


Putting together the expressions for F + and F , we finally end up with the vacuum 
eigenvalue of the transfer matrix when 0 < Ima < 7 t/2: 


A(a) = exp[F + (a)] + exp[F (a)] = 


= 2 


r + uj — £7re _2£ “ 
r (| + uj + £7re _2e “) 

2 7Te 2iriP(l)e~ 


e 2 ^(i) e 2ea cosh(i7r 2 £e 


2 _ 2 sol 


— emu = 


r (| + w + Tre 2ea ) r^-u + ne 2e ") 


(5.93) 
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This gives the expression for the vacuum eigenvalue of our conformal transfer matrix 
at /3 2 = 1/2 and in the region 1 (15.911 of the spectral parameter a; it is, of course, an 
entire function of A 2 = e 2 “. 

In region 2 (ODD , he. —7 t/ 2 < Ima < 0, the calculations are carried on in a similar 
way. We only remark that the computation of F + (a) is made starting from formula 
(Oil and that, because of the relation (IQl 


F (a) = — F + (a + ij) = —F + ^a + + iivr^J , 


(5.94) 


the evaluation of F (a) comes from formula on . The final result is: 


F ± (a) = 


-e in T ( - + uj + ene 


— 2ea 


+ e in r 


+ u — Eire 2ea ) + 


+ 27 T 0 (l)e 2£a ±m 2 ee 2ea =F emu . 


2 _ 2 ea 


(5.95) 


These expressions coincide in form with (15.881 15.921) . therefore expression (15.881) for 
F + (a ) and expression (I5.92|) for F~(a) are valid in all the strip —7 t/ 2 < Imcc < 7 t/ 2. 
Since F ± (a) are periodic with period m Oil) and since (PKi on show the same 
periodicity, we conclude that in all the complex plane the expressions of F + {a ) and 
F~(a) at /3 2 = 1/2 are given respectively by (15.8811 and (15.9211 . It follows that the 
vacuum eigenvalue of the transfer matrix in all the complex plane at (3 2 = 1/2 is given 
by expression on . 

In the case e = — 1 expression on becomes 


A(a) 


2 7re 2 7 n/’( 1 )e 2 “ 

T (i +cn + 7re 2 “j T — uj + 7re 2 ") 


(5.96) 


which means that, after the usual identifications A 2 = e 2 “ and (Remark 4) a; = 2p, 
A(a) coincides with the expression for the transfer matrix at j3 2 = 1/2 given in formula 
(4.32) of [2]. 

Expression (E2SD satisfies a peculiar functional equation, which identifies as its 
solutions all the transfer matrix eigenvalues (at the free fermion point), once we assume 
the solution to be entire. This functional equation is nothing but the renormalised 
version of the T-Q relation which holds for 0 < p 2 < 1/2 PUG]. Indeed, we can 
decompose A ± (a) in the following form 

exp[F ± (a)] = exp(±i7r 2 ee _2£ “ =F sinco )— ^ 2 , (5.97) 

q/(a) 

where we have defined 


Q(a) 


T (A + u + 7ree- 2£Q ) 
28 


exp [— 2sa ] . 


(5.98) 























Since 


Q [a 


17V 

2 


Q[ 


(X + 


Q(a) Q(a ) 

expression (15.931) can be written as follows 


( X7V \ / \ 

a + — J + exp(— iTT 2 ee~ 2£a + £incu)Q (a —— J . 

(5.99) 

Therefore, assuming that our model at £ = —1 is a lattice regularisation of the contin¬ 
uous theory in 2], we have proved the functional equation (4.27) of [2]. This equation 
has been proved there as exact at first order in perturbation theory and then used to 
derive (15.961) . 


Finally, we present the calculation, at the free fermion point, of the local integrals 
of motion. We can specialize formula flElBD - with A chosen as in Remark 4 - or, alter¬ 
natively, expand asymptotically the logarithm of (15.931) . taking into account relations 

(IQK1ICT1ICT1) : 


9 F 00 

lnA(a) = [27t(1 + -0(1) — ln7r) + Ansa] e~ 2sa -- e 2ea(2n+1 )c n+1 I 2n +i , (5.100) 

* n=0 


where 


and 


Cn+1 


hn+l — —£ 


2 n 1 
7r 2n (n + l)(2n + 1) 

^+T- B2n+2 (2 + W ) ’ 


(5.101) 

(5.102) 


where B 2m (x), m = 1,2,..., are the Bernoulli polynomials |26j . The previous formula 
(15.1021) has been already found using the free fermion basis in [521J . 


6 Summary and conclusions 

In the whole interval 0 < (3 2 < 1, we have determined the nonlinear integral equation 
for the a)-vacuum states of the left and right (m)KdV models, whose discretisations 
are dynamically twisted spin —1/2 XXZ chains. The mapping /3 2 —> 1 — (3 2 connects 
the (m)KdV nonlinear integral equation with the nonlinear integral equation coming 
from the twisted spin +1/2 XXZ chain. We have identified conformal primary states 
in those vacua and the treatment of excitations - as well as off-critical behaviour - 
can be derived along the lines of H]- We have given an expression for the eigenvalues 
of the transfer matrix in terms of solutions of the nonlinear integral equation in the 
whole interval 0 < /3 2 < 1, and we have found a rich analytic structure, which in 
particular implies different asymptotic expansions. From the latter the expressions 
of the local integrals of motion shared with KdV theory have been extracted. At 
the free fermion point the aforementioned eigenvalues are evaluated exactly, again in 
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agreement with 0. Therefore, we have been naturally led to conclude that our model 
furnishes an effective framework for calculating interesting quantities of both quantum 
KdV theory and Conformal Field Theories. Moreover, the exact values of the transfer 
matrix eigenvalues may be useful in VT - symmetric Quantum Mechanics applications, 
for they have been proved to be spectral determinants in In this respect, the 

exact solution of the theory given in the limit of infinite twist (u —> +oo) may be of 
particular interest in many applications and developments. Finally, we would like to 
highlight that the second paper will concern the off-critical theory approach proposed 
in [BJ within this framework, showing the generality of the methodology. 
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